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FOREWORD 


A  computer  program  is  presented  which  is  capable  of  modelling  the  geometry 
and  performance  of  the  water  entry  cavity  produced  by  the  entry  of  axially 
symmetric  projectiles.  The  program  is  based  on  an  existing  theoretical  cavity 
flow  model  modified  to  account  for  atmospheric  surface  pressure.  The  computer 
model  is  shown  to  satisfactorily  predict  cavities  produced  by  vertical  and 
oblique  entries  (including  high-speed  entries)  from  normal  air  of  right  circular 
cylinders  and  truncated  ogives. 
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INTRODUCTION 


The  iirst  part  ot  this  report  presents  the  final  version  of  a  computer 
program  which  is  capable  of  modelling  the  behavior  of  Che  water  entry  cavity 
produced  at  vertical  water  entry.  The.  program  is  baaed  on  a  hydraulic  cavi  ty 
flow  model  proposed  some  years  ago  by  G.  Birkhoff  and  K.  Isaacs^  which  has 
been  modified  to  include  the  effects  of  an  atmospheric  surface  pressure. 

The  second  part  of  this  report  presents  a  summary  <  *'  the  results  of  a 
theoretical  and  experimental  study  of  the  water  entry  cavity.  In  this  study, 
observed  cavities  are  compared  with  cavities  predicted  by  the  computer  program 
described  in  Part  I.  Various  water  entry  conditions,  covering  a  wide  range  of 
entry  velocities,  were  investigated. 


Ifiirkhoff,  G.  and  Isaacs,  R*,  "Transient  Cavities  in  Air~Water  Entry," 
NAVORD  1490  (1951). 
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PART  1  -  CAVITY  MODELLING  CODE 


liACKGROUWD 


Birkhoff  and  Isaacs  developed  two  flow  models  ~  one  which  applies  Co 
vercicfll  entry  and  the  other  to  oblique  entry.  The  discussion  here  will  be 
limited  to  the  vertical  entry  model  on  which  the  program  is  based.  The  original 
model  is  applicable  to  the  prediction  of  the  open  cavity  produced  by  the  water 
entry  of  axially  aymmetric  missiles  from  either  a  vacuum  or  greatly  reduced 
pressure  atmosphere  (Abelson^  found  that  it  also  gives  good  shape  predictions 
of  cavities  produced  by  normal  air  entry  up  until  the  time  of  pullaway). 

Birkhoff's  vertical  entry  model  is  based  on  the  hypothesis  that  flow 
"streamlines  lie  on  concentric  spherical  surfaces  centered  at  the  point  of 
impact."  Accordingly,  the  liquid  medium  can  be  thought  of  as  being  subdivided 
into  infinitesimally  thin  concentric  spherical  shells  or  "laminas"  centered  at 
the  entry  point,  as  shown  in  Figure  1.  As  the  missile  passes  a  lamina  the  fluid 
within  is  displaced,  or  forced,  radially  outward  as  shown  in  Figure  2.  The 
inside  edges  of  the  open  laminas  together  form  the  cavity  wall  outline.  The 
shells  derive  their  initial  kinetic  energy  from  the  loss  of  kinetic  energy  of 
the  passing  missile,  resulting  from  the  cavity  drag.  The  cavity  drag 
coefficient  of  the  missile  is  assumed  to  remain  constant  resulting  in  an 
exponential  velocity  decay  with  depth 


V 

P 


r 

P 


U) 


where  ■ 

Vp^  *=  missile  entry  velocity 
rp  ••  missile  depth 
a  ■  retardation  coefficient 


2  Abel  son,  h. ,  "The  behavior  of  the  Cavity  Formed  by  a  Projectile  Entering  the 
Water  Vertically,"  Ph.D.  Thesis,  Mechanical  Engineering,  Univ.  of  Md.,  1969. 
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Considering  a  typical  lamina  located  a  distance  r  Irom  the  entry  point) 
Figure  3,  only  the  location  of  the  inside  edge  of  the  displaced  shell,  as 
defined  by  r  and  the  angle  9^,  is  required  in  order  to  determine  the  position 
of  the  cavity  boundary.  By  treating  the  shell  as  a  kind  of  "fluid  pendulum"  and 
utilizing  the  principle  of  conservation  of  mechanical  energy,  the  differential 
equation  relating  the  angle  8^  to  time  can  be  shovm  to  be^>- 


CD 

1 

o 

• 

r 

•s^ 

-  V’ 

r  J  sin^  0^  In  T 

(2  - 

cos  0a) 

(1  +  cos  0a) 

cos 

0a  ^  “ 

cos  9a) 

(2) 


where  the  expression  for  cos  0^  ia*^ 


cos 


1 

2 


0.  Op  Vp2  e'Zar 

TT  r3  p„  g 


(3) 


The  angle  9^  is  the  maximum  value  of  the  displacement  6^  corresponding  to 
the  peak  of  the  lamina's  swing*  The  initial  condition  can  be  found  by 


R 

(o)  “  arctan  -E 

**0  ®  r 


U) 


where  0^  .  the  initial  shell  angle,  is  tsken  by  Birkhoff^  to  be  the  angle 
subtended  at  Che  point  of  entry  by  the  zone  of  separation  of  the  water  from  the 
missile  surface.  This  would  imply  Chat  Kp  should  be  taken  as  the  radius  of 
the  missile  at  the  zone  of  separation*  For  certain  missile  configurations  it  is 
obvious  where  the  separation  occurs  and  thus  what  Kp  should  be  (for  example, 
for  a  right  circular  cylinder  only  the  nose  flat  is  wetted  hence  Rp  would  be 
taken  as  half  the  nose  flat  diameter)*  For  certain  other  missile  configurations, 
on  Che  other  hand,  the  location  of  separation  n\ay  not  be  obvious  and  may  be 
difficult  to  determine*  In  such  cases  it  may  be  necessary  to  estimate  Rp  (for 
example  in  Abelson's  cavity  study^  which  involves  sphere  entries,  Rp  is 
taken  as  Che  transverse  radius  of  the  missile  body  at  its  center  of  gravity, 
i*Q<,  halt  the  diameter  for  a  sphere)  •* 


^See  footnote  1  on  page  7. 

^Abel  son,  H. ,  "Cavity  Shapes  at  Vertical  Water  Entry  -  A  Comparison  of 
Calculated  and  Observed  Shapes,"  NOLTR  67-31  (1967). 

*In  reference  (3)  the  expression  for  cos  0[^  is  in  error  in  that  the  first  term 
in  the  brackets  in  Equation  (3)  appears  as  a  "l"  instead  of  a  "2"  as  it  should 
be.  The  author,  apparently  through  an  oversight,  failed  to  account  for  the  "at 
rest"  potential  energy  of  the  shell  which,  because  of  the  reference  datum 
selected,  is  not  equal  to  zero.  If  this  at  rest  term  is  included  in  tlie 
reference  (3)  derivation  the  correct  expression  for  cos  with  the  "2"  in  the 
brackets  Instead  of  the  "1",  is  obtained. 
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By  differentiating  (2)  with  respect  to  time  the  expression  for  the  angular 
acceleration  is  obtained: 


0  - 

a  r 


1  +  (cos  6  -  cos  6  )  Z 

_ a _ m 

sin  0  In  (k} 

La  J 


where : 


2-4  cos  0, 


2  cos  6 


COB  G  (cos-^  0  -  2  coB^  0  -  cos  9  +  2)  In  (jO  sin^  0 

a  a  a  a  a 

(2  -  cos  0.)  (1  +  cos  0„) 

at  •  I  ^  ^ 

cos  0a  (1  -  cos  8a^ 

MODIFICATIONS  TO  ORIGINAL  FLOW  MODEL. 

The  foregoing  equations,  which  describe  Birkhoff’s  original  model,  do  not 
account  for  pressure  either  in  the  cavity  or  at  the  water  surface  since  they 
were  assumed  to  be  zero.  However,  as  mentioned  earlier,  the  model  as  used  in 
the  program  has  been  altered  to  include  the  effects  of  surface  pressure  since 
this  modification  is  fairly  straightforward,  providing  the  surface  pressure  is 
regarded  as  constant.  Modifications,  on  the  other  hand,  to  include  the  effects 
of  pressure  in  the  cavity  on  cavity  behavior  are  not  a  simple  matter  at  all.^ 
Thus,  in  the  program  cavity  pressure  is  assumed  to  equal  zero. 


When  the  effect  of  surface  pressure,  Pq,  is  included  it  is  shown  in 
Appendix  2  that  the  governing  differential  equation  for  the  displacement  angle 


0g  becomes; 


^a  -  i 


4 

rg  +  (Po/P„)_ 

COtt  0  “  COfl  0 

a  ta 

2  ^.2  Q  l-vl 

where,  as  in  the  original  model,  the  initial  condition  is 


0fl  "  Cretan 
“  0 


with  corresponding  relations  for  Gg  and  cos  0,^  as 

-2  frg  +  (P„/p  )1  fl  +  0.  -  COB  9  )  (Z)' 

fl  - _ I - 1-JLl _ J. - - 


sin  6  In 
a 


^See  footnote  2  on  page  9. 
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and 


cos  6- 


am 

2  - 

TT  r2  (r 


P  ® 
0 _ 


-2ar 


Pw  B  +  Po) 


(8) 


It  can  be  seen  t'\ut  by  setting  Pq  equal  to  zero,  Equations  (6),  (7)  and 
(8)  reduce,  to  Equations  (2),  (3)  and  (5). 


Clearly  the  differential  equations  for  the  angle  0^ 
and  the  modified  model  are  non-linear  and  must  be  solved 
using  numerical  integration  since  analytic  solutiona  for 
yet  to  be  found.  The  program,  in  fact,  solves  Equations 
numerical  integration. 


for  both  the  original 
for  0^  versus  time 
these  equations  have 
(6) ,  (71  and  (8)  using 


CAVITY  MODELLING 


It  should  be  remembered  that  these  lamina  equations  apply  to  a  single  lamina 
at  a  distance  r  from  the  entry  point.  Hence,  the  motion  of  a  liinina  at  any  depth 
can  be  found  independent  of  the  laminas  at  the  other  depths.  For  this  reason  the 
entire  cavity  may  be  modelled  satisfactorily  with  the  use  of  a  relatively  few 
number  of  laminas  located  at  selected  depths,  as  shown  in  Figure  4.  The  cavity 
outline  at  any  given  time  may  be  obtained  easily  by  plotting  the  positions  of 
the  edges  of  the  laminas,  then  connecting  the  appropriate  points  to  get  the 
cavity  outline. 


Another  less  accurate,  but  more  convenient  method  of  describing  cavity 
shape,  which  will  represent  the  cavity  for  all  times,  involves  plotting  the 
cavity  diameter  at  the  selected  depths  versus  time  after  entry.  The  result  is  a 
set  of  arc-shaped  curves  like  the  ones  shown  in  Figure  8.  Obviously,  from  the 
figure,  deeper  laminas  start  their  motion  at  later  times.  It  can  easily  be 
shown  that  when  cavity  drag  coefficient  is  constant  the  time  for  the  missile  to 
reach  a  lamina  at  depth  r  is: 


t  =  (e°“^  -  l)/avp 


From  Figure  5,  the  diameter  of  the  circular  opening  in  the  lamina  at  a  given 
time  is  found  by  the  relation 

D  “  2r  sin  0^  (9) 

There  is  an  approximation  being  made  here  however.  Although  the  diameter 
computed  from  Equation  (9)  is  exact,  it  is  actually  the  diameter  of  the  cavity 
not  at  the  depth  r  but  at  the  more  shallow  depth  r(l  -  cos  Sg)  since,  from 
Figure  5,  the  edge  of  the  lamina  rises  up  as  the  laminn  swings  out.  For  small 
enough  values  of  the  angle  0g  this  should  present  no  problem.  However,  even 
for  fairly  large  values  of  0g,  the  approximation  should  still  be  reasonable 
since,  in  general,  the  cavity  wall  is  nearly  vertical  meaning  the  change  in 
diameter  with  depth,  at  a  given  time,  is  not  that  great.  A  final  consequence  of 
using  this  technique  is  that  at  a  given  depth  r,  the  cavity  diameter  will  never 
exceed  the  limiting  value  2  r. 
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The  program  is  designed  to  make  it  convenient  to  implement  either  of  the 
two  modelling  techniques  mentioned.  The.  output  data  are  in  the  form  of 
"Spherical  Lamina  Performance  Data"  tables  where  each  table  corresponds  to  a 
lamina  at  some  user  specified  depth  and  contains  information  concerning  depth  r, 
0a,  and  Cavity  Diameter  versus  absolute  time  after  entry  as  well  as  other 
auxilliary  information  such  as  current  model  depths  and  angular  velocities  and 
accelerations  associated  with  the  shell  displacement  angle  6g. 

MINIMUM  DEPTH  RESTRICTIONS.  The  selection  of  depths  at  which  to  model  the 
cavity  cannot  be  made  entirely  arbitrarily  since  the  flow  model  equations  break 
down  in  the  region  near  the  water  surface. 


For  any  given  entry  there  is  a  lamina,  located  at  a  critical  depth, 
designated  r^j.^,  which  acquires  just  enough  energy  from  the  passing  missile  to 
rise  up  to  where  is  equal  to  90°  (i.e.,  entry  conditions  are  such  that 
for  this  particular  lamina  6^  is  equal  to  90°).  Needless  to  say,  all 
laminas  at  depths  less  than  r^.^.  acquire  sufficient  energy  to  swing  up  beyond 
the  point  0g  ■■  90°.  In  physical  terras  such  a  lamina  will  continue  to  swing 
upward  until  it  flows  into  itself  at  the  point  0a"  90°  as  depicted  in 
Figure  6.  By  setting  cos  9 m  in  Equation  (6)  equal  to  zero  (corresponding  t'. 

Qp,  “  90°)  and  replacing  r  with  rcj.  ,  the  following  formula  is  obtained 
which  can  be  solved  by  tr ial-and-error  for  the  first  critical  depth  ; 


0  = 


1 

2 


2 


01  m 


2  -2arr 


IT 


1 

8  * 


Inspection  of  Equations  (6)  through  (8)  will  show  that  it  is  safe  to  model 
the  cavity  at  depths  greater  than  r^ri,  which  assures  that  the  displacement 
angle  Bg  will  not  exceed  90°.  Attempting  to  model  at  depths  equal  to  or  less 
than  rcri  will  result  in  a  singularity  in  Equations  (6)  through  (8)  when  0^ 
reaches  90°,  although  modelling  up  until  this  condition  occurs  may  be  possible. 


A  second- critical  depth,  designated  ,  occurs  nearer  to  the  surface 

than  and  is  the  depth  below  which  Equation  (8),  for  cos  0,^,,  breaks 

down  completely.  Attempting  to  model  at  depths  less  than  will  result  in  a 
value  for  cos  0^  which  is  less  than  -1.  Even  so,  as  in  the  case  of  the  first 
critical  depth.  Equations  (6)  through  (8)  still  "work"  until  0g  reaches  90° 
(where  the  singularity  occurs)  although  it  is  not  known  at  this  time  if  the 
results  produceil  will  be  meaningful. 


Figure  7  summarizes  the  minimum  depth  restrictions.  Since  these 
restrictions  apply  to  depths  in  a  small  region  near  the  water  surface  they 
should  not,  in  general,  interfere  with  the  modelling  of  the  overall  cavity. 


PROGRAM  DESCRIPTION 


The  computer  program  is  written  in  standard  FORTRAN 
appears  in  Appendix  1.  The  first  part  of  the  listing  is 
manual."  To  demonstrate  the  use  of  the  program  a  sample 
carried  out  in  the  next  article.  A  brief  description  of 
here . 


and  a  listing  of  it 
essentially  a  "user's 
problem  has  been 
the  program  is  given 
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Basically  the  program  solves  the  lamina  motion  Equation  (6),  (7)  and  (8) 
for  3^  with  time  for  the  lamina  at  each  user  specified  depth  using  the  three 
term  taylor  series: 


0 


a(i  +  1) 


°“(i)  *  *•(!) 


Currently,  the  maximum  allowable  number  of  depths  which  may  be  specified  is 
10,  however  this  number  may  be  increased  simply  by  increasing  the  Fortran 
dimension  parameters  for  the  variables  DEPTH,  DELTA  and  NSTORE  as  required.  The 
end  result  of  the  computations  for  a  lamina  at  a  particular  depth  is  the  single 
"Spherical  Lamina  Performance  Data"  table  described  in  the  previous  section. 
There  are  no  "Max  Time"  or  "Max  Iterations"  parameters  to  bother  with  since, 
when  the  motion  of  a  particular  lamina  is  completed,  the  sign  of  the  angle  6^ 
goes  negative,  and  the  computations  for  that  lamina  are  terminated 
automatically.  A  sample  of  the  data  tables  appears  in  Appendix  1  following  the 
program  listing.  These  tables  are  the  results  for  the  sample  problem  rlescribed 
next . 


SAMILE  PROBLEM 


The  program  was  used  to  model  the  cavity  produced  by  the  water  entry  of  a 
steel  right  circular  cylinder,  1.50"  diameter  by  2.15"  long,  from  normal  air,  at 
an  entry  speed  of  250  feet/second.  The  cavity  was  modelled  with  laminas  at 
depths  of  10,  20  and  30  calibers.  Several  test  runs  were  carried  out,  using 
progressively  smaller  values  for  the  time  step,  until  convergence  to  the 
solutions  for  all  three  laminas  was  obtained.  The  results  of  the  final  run  can 
be  found  in  Appendix  1  following  the  program  listing.  The  first  page  of  the 
output  is  a  Reference  Information  table  listing  data  describing  the  entry 
conditions.  Following  the  Reference  Information  are  the  three  "Spherical  Lamina 
Performance  Data"  tables  for  the  laminas  at  each  of  the  selected  depths.  The 
data  in  these  tables  were  used  to  generate  the  curves,  in  Figure  8,  of  cavity 
diameter  history  at  each  of  the  three  depths. 
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PART  2  -  CAVITY  STUDY 


BACKGROUND 

G.  Birkhoff^  compared  cavities  predicted  by  liis  original  model  tor 
vertical  entry  with  experimental  data  concerning  a  one  inch  diameter  steel 
sphere  entering  water  vertically  from  a  reduced  pressure  atmosphere  at  a 
velocity  of  73  feet/second  with  a  surface  pressure  of  1/27  atmosphere* 

Abelson^  compared  the  cavities  predicted  by  Birkhotf's  model  to  those 
generated  by  the  vertical  entry  of  spheres  into  water  from  reduced  pressure 
atmospheres.  The  spheres  had  diameters  of  either  0.5  or  1.0  inches  and  entry 
speeds  of  between  62  and  69  feet/sccond. 

Initially,  Birkhoff’s  original  model  was  applied  only  to  open  cavities 
produced  by  entry  from  a  vacuum  or  reduced  pressure  atmosphere.  Then 
Abelson^  compared  the  model  to  the  cavity  produced  by  a  normal  air  entry  and 
found  that  theoretical  and  observed  cavity  shapes  agreed  well  up  until  the  time 
of  pullaway.  Abelson  then  modified  Che  model  to  include  the  effects  of  both 
constant  surface  pressure  and  a  measured  cavity  pressure.  He  Chen  compared  the 
cavities  predicted  by  the  altered  model  to  the  observed  normal  air  entry 
cavities  and  found  that  the  modified  model  predicted  cavity  shapes  which  were, 
in  fact,  considerably  smaller  in  size  than  both  the  observed  cavities  or  those 
predicted  by  the  original  model.  The  fact  that  Birkhoff's  original  model, 
unaltered  for  system  pressure,  satisfactorily  predicts  cavities  produced  by  both 
normail  and  reduced  pressure  entries  suggested  to  Abelson  that  cavity  shape  might 
be  practically  independent  of  surface  or  cavity  pressure,  at  least  up  to  the 
time  of  pullaway.  In  an  attempt  to  confirm  this  hypothesis,  Abelson  compared 
the  cavities  produced  by  sphere  entries  from  both  normal  and  reduced  pressure 
atmospheres  and  found  that  "the  cavities  were  identical  until  surface  closure 
and  then  began  to  deviate  increasingly  until  pullaway,  with  the  cavities  of  the 
reduced  pressure  entries  slightly  larger."  As  a  result  of  his  observations, 
Abelson  concluded  that  the  cavity  shape  is  in  fact  practically  independent  of 
surface  or  cavity  pressure  up  to  the  time  of  pullaway  and  hence,  that  Birkhoff's 
original  model  is  the  one  to  use  to  satisfactorily  predict  cavity  shape  and  that 
"any  modification  of  the  model  to  include  system  pressure  invalidates  it." 


^See  footnote  1  on  page  7. 

a 

See  footnote  2  on  page  9. 

2 

See  footnote  3  on  page  10, 
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EXPERIMENTAL/THEORETICAL  STUDY 

Comparison!  of  Birkhoff's  vertical  entry  model  to  experiment  are  limited 
mostly  to  shots  involving  spheres  entering  at  relatively  low  velocities  (below 
73  feet/second).  An  exception  is  Abelson's  comparison  of  the  model  to  the  case 
of  a  water  entry  involving  a  3  inch  diameter  cylindrical  missile  having  a  140 
degree  conical  nose,  entering  at  200  feet/aecond  from  normal  air. 2  Jn  view  of 
these  limitations  it  was  felt  that  it  would  be  useful  to  compare  both  the 
original  model  and  the  version  modified  for  surface  pressure  (in  spite  of 
Abelson's  conclusion  that  inclusion  of  system  pressure  invalidates  the  model)  to 
entries  involving  higher  velocities  and  model  configurations  other  than  spheres. 

Specifically,  theoretical  cavities  predicted  by  the  models  were  compared  to 
observed  cavities  from  a  number  of  shots  involving  normal  air  pressure  entries 
of  right  circular  cylinders  of  various  diameters  and  aspect  ratios  and  also  to 
shots  involving  cylindrical  missiles  of  various  sizes  having  truncated  ogive 
nose  configurations. 

One  source  of  experimental  data  was  photographic  records  from  vertical 
water  entry  testa  conducted  in  the  NSWC  Pilot  Tank  Facility  in  1967  under  the 
supervision  of  Dr.  A.  May.^  The  shots  chosen  for  this  study  involved  right 
circular  cylinders  and  stepped  right  cylinders.'*  Model  and  test  condition  data 
for  the  shots  chosen  from  May's  test  for  this  study  are  given  in  Table  1.  Nose 
flat  diameters  ranged  from  0.9  to  1.5  inches  and  entry  velocities  ranged  from 
175  to  415  feet/second. 

High-speed  motion  pictures  from  another  water  entry  test  program,  conducted 
in  Che  NSWC  Undersea  Weapons  Tank  around  1959,  provided  a  second  source  of 
experimental  daCa.^  This  test  scries  involved  high  speed  oblique  entries  of 
truncated  ogives  and  right  circular  cylinders.  Although  these  shots  involved 
oblique  entries  it  was  felt  that,  for  the  depths  of  interest  in  this  study 
(below  50  calibers**),  the  entry  angles  of  45,  47  and  70  degrees  (from 
horizontal)  were  steep  enough  that  changes  in  the  cavity  geometry,  as  compared 
to  a  vertical  entry,  would  not  be  appreciable.  For  the  shots  chosen  from  the 
series  for  this  study  the  entry  model  nose  flat  diameters  ranged  from  0.707  to 
1.570  inches  and  entry  velocities  ranged  from  426  to  1175  feet/second.  Model 
and  test  condition  data  are  given  in  Table  2. 


^See  footnote  2  on  page  9. 

^May,  A.,  "The  Cavity  After  Vertical  Water  Entry,"  NOLTR  68-114  (1968). 

^May,  A.  and  Hoover,  W.  R. ,  "A  Study. of  the  Water-Entry  Cavity," 

NOLTR  63-264  (1963). 

*For  the  two  shots  involving  the  stepped  cylinders,  the  stepped  shoulder 
remained  within  the  cavity  and  was  not  observed  to  interfere  in  any  way 
with  cavity  production. 

**Ba8ed  on  nose  flat  diameters. 
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TABLE  1 


MODEL  AND 

TEST  CONDITION  DATA: 

PILOT  TANK 

SERIES 

SHOT  # 

MODEL 
CONFIG . 

NOSE  FLAT 

DIA. 

MODEL 

LENGTH 

MODEL 

WT. 

MAT'L. 

(in)  (in)  (lbs) 


1419 

RIGHT  CIRC.  CVL. 

1.000 

8.25 

1.834 

STEEL 

1711 

RIGHT  CIRC.  CVL. 

1.500 

2.00 

2.143 

UNKNOWN 

1438 

RIGHT  CIRC.  CYL. 

1.500 

4.39 

2.195 

STEEL 

1434 

RIGHT  CIRC.  CYL. 

1.500 

2.15 

1.075 

STEEL 

1474 

SINGLE  STEPPED  CYL.* 

0.900 

10.10** 

1.818 

STEEL 

1476 

SINGLE  STEPPED  CYL.* 

0.900 

10.10** 

1.818 

STEEL 

ENTRY 

DRAG 

CAVITY 

ENTRY 

SHOT  # 

VELOC. 

AREA 

DRAG  COEFF. 

ANGLE 

( ft/ sec)  (sq  ini) 


1419 

175.0 

.785 

.807 

VERTICAL 

1711 

200.3 

1.767 

.807 

VERTICAL 

1438 

245.0 

1.767 

.807 

VERTICAL 

1434 

250.0 

1.767 

.807 

VERTICAL 

1474 

320.0 

0.636 

.807 

VERTICAL 

1476 

415.0 

0.636 

.807 

VERTICAL 

*Stepped  cylinder  dimensions:  0.90"  dia  x  0.975"  long  nose  cylinder; 

1.50"  dia.  X  3.308"  long  body  cylinder. 

**Length  for  an  equivalent  steel  unstepped  cylinder  (i.e.,  one  having  the 
same  weight  and  nose  flat  diameter  of  the  stepped  cylinder). 
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TABLE  2 


MODEL  AND  TEST 

CONDITION 

DATA:  UNDERSEA 

WEAPONS  TANK 

SERIES 

MODEL 

NOSE  FLAT 

MODEL 

MOD 

SHOT 

1 

CONFIG . 

DIA.* 

LENGTH 

WT. 

MAT'L. 

(in) 

(in) 

(lbs) 

18A 

TRUNCATED  OGIVE 

0.707" 

NOT  AVAILABLE 

1.599 

NOT  AVAILABLE 

192 

TRUNCATED  OGIVE 

0.864 

NOT  AVAILABLE 

1.600 

NOT  AVAILABLE 

205 

TRUNCATED  OGIVE 

.774 

NOT  AVAILABLE 

1.669 

NOT  AVAILABLE 

193 

TRUNCATED  OGIVE 

.864 

NOT  AVAILABLE 

1.600 

NOT  AVAILABLE 

63 

TRUNCATED  OGIVE 

.785 

NOT  AVAILABLE 

1.652 

NOT  AVAILABLE 

172 

RIGHT  CIRC.  CYL. 

1.570 

NOT  AVAILABLE 

1.537 

NOT  AVAILABLE 

ENTRY 

DRAG 

CAVITY 

ENTRY 

SHOT 

I 

VELOC. 

AREA 

DRAG  COEFF. 

ANGLE** 

, 

(f t/aec) 

(sq  in) 

(deg) 

18'+ 

426.0 

.393 

.800 

47.0 

192 

548.0 

.586 

.800 

47.0 

205 

599.0 

.471 

.800 

47.0 

193 

632.0 

.586 

.800 

47.0 

63 

946.0 

.484 

.800 

70.0 

172 

1175.0 

1.936 

.800 

45.0 

*For  these  models  full  body  diameters  were  equal  to  1.57  inches. 
'*'*Angle  between  theoretical  gun  line  and  water  surface. 
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SUMMARY  OF  RESULTS 

All  together)  twelve  water  entries  were  modelled  with  results  presented  in 
the  form  of  cavity  diameter  histories  shown  in  Figures  9  through  23.  The  curves 
in  Figures  9  through  20  compare  observed  cavities  with  cavities  predicted  by  the 
altered  theoretical  model  which  includes  effects  of  surface  pressure.  The 
remaining  curves  in  Figures  21  through  23  compare  observed  cavities  with 
cavities  predicted  by  Birkhoff s  original  model.  Each  graph  contains  cavity 
diameter  time  histories  for  both  predicted  and  observed  cavities  at  selected 
depths.  The  data  in  nearly  all  the  plots  extend  well  beyond  the  time  of 
pullaway  (the  only  exception  is  shot  iH72  where  a  deep  closure  occurred  and  no 
pullaway  was  observed).  Time  data  were  obtained  from  timing  marks  present  on 
th  e  film. 

For  the  curves  in  Figures  9  through  23,  corresponding  to  the  study  with 
surface  pressure,  the  fallowing  observations  were  made: 

1.  In  general  the  model  altered  for  surface  pressure  appears  to 
underestimate  the  cavity  diameter,  hence  cavity  size.  This  observation  should 
be  expected  to  some  degree  since  cavity  presure  is  assumed  to  be  zero.  The 
exceptions  are  shots  #63  and  #1476.  In  both  these  cases  curves  for  observed  and 
predicted  cavities  agree  closely  for  all  depths  for  roughly  the  first  half  of 
the  curves,  corresponding  to  cavity  growth,  then  for  the  last  half, 
corresponding  to  cavity  collapse,  the  model  overestimates  the  diameters. 

2.  Model  correctly  predicts  a  deep  closure  type  of  collapse  (collapse 
from  Che  bottom  up)  for  shot  #172  and  base  closure  type  collapse  (collapse  from 
top  of  the  cavity  downward)  for  the  remaining  shots. 

3.  Theoretical  "lag  times"  (times  needed  for  missile  to  reach  selected 
depths)  which  have  nothing  Co  do  with  the  model  but  result  from  Che  assumption 
of  constant  cavity  drag  coefficient,  agree  well  with  observed  lag  times. 

Following  the  cavity  diameter  histories  is  a  set  of  curves  in  Figure  24 
showing  the  results  of  a  parametric  study  of  Che  surface  pressure  Pq  to 
determine  Che  sensitivity  of  cavity  size,  at  a  representative  depth  of  20 
calibers,  to  variations  in  surface  pressure. 

The  cavity  diameter  histories  in  Figure  25  are  for  a  theoretical  cavity 
predicted  by  Birkhoff's  original  model.  It  is  included  to  show  how  the  typical 
theoretical  cavity  behaves  since,  for  the  study  involving  the  original  model  in 
Figures  21  through  23  only  a  small  portion  of  the  theoretical  curves  could  be 
shown  on  the  graphs  because  of  their  extremely  long  periods,  as  compared  to  the 
experimental  curves.  The  diameter  histories  in  Figure  25  (essentially  an 
extension  of  the  theoretical  curves  in  Figure  15  for  shot  #1474)  indicate  that  a 
deep  closure  occurs  at  a  depth  of  140  calibers. 

RECOMMENDATIONS 

Ihe  following  recommendations  suggest  areas  to  consider  for  future  work  and 
items  which  may  be  of  concern  in  a  more  involved  study: 
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1.  Some  years  ago,  according  to  Hi  Ki  Steves,  there  were  ahota  done  in 
the  NSWC  Hydro-Balliatics  Tank  Facility  involving  high  speed  vertical  entries 
which  are  recorded  on  filmi  It  may  be  useful  to  see  how  well  cavities  produced 
by  the  model  compare  with  these  cavities  since  the  high  speed  shots  considered 
in  this  study  involved  oblique  entries  whereas  the  flow  model  is  intended  to  be 
applied  to  vertical  entries. 

2.  Reread  the  photographic  records  for  the  vertical  entries  in  Figures 
21  through  23  to  determine  approximately  when  the  time  of  pullaway  occurs  and 
indicate  these  times  on  the  graphs.  Then  compare  the  results  in  these  curves 
with  Abelson's  suggestion  that  cavity  shape  can  be  satisfactorily  predicted  up 
until  pullaway  for  normal  air  entries  using  Birkhoff's  original  model. 

3.  Modif icati  .1 1'  nould  ba  made  to  the  current  version  of  the  cavity 
modelling  code  to  enable  to  handle  the  problems  associated  with  modelling  the 
cavity  in  the  critical  mlni.mum  depth  region.  Currently  no  such  capability 
exists  in  the  code,  consequently  programming  difficulties  related  to  the 
numerical  complications  described  earlier  will  arise  if  an  attempt  is  made  to 
model  in  this  region.  At  the  very  least,  changes  could  be  made  to  the  code  so 
that  it  will  simply  ignore  specified  depths  which  fall  within  the  critical 
region  and,  in  addition,  issue  a  message  warning  the  user  of  the  condition. 

More  sophisticated  modifications  could  allow  some  limited  modelling  within  this 
critical  region. 


20 


NSWC  TR  81-69 


LINE 


FIGURE  1  LIQUID  MEDIUM  SHOWN  SUB-DIVIDED  INTO 
CONCENTRIC  SPHERICAL  SHELL  OR 
"LAMINA"  REGIONS 


l.q.f. 


FIGURE  2  DISPLACEMENT  OF  FLUID  SHELLS  TO  F  ORM  CAVITY 
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FIGURE  4  MODELLING  THE  CAVITY  WITH  THE  USE  OF  LAMINAS  AT 
SELECTED  DEPTHS 


FIGURE  B  CAVITY  DIAMETER  APPROXIMATION.  DIAMETER  AT  DEPTH  r 
APPROXIMATED  BY  PROJECTING  D  DOWN  TO  DEPTH  r 


FIGURE  6  FLUID  SHELL  IMPACTING  ITSELF 
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L.O.F. 


FIGURE  7  SUMMARY  OF  CRITICAL  MINIMUM  DEPTH  RESTRICTIONS 
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SAMPLE  RUN 

ENTVEL-  260  Ipi 
QUN  ANG  -  VERTICAL 

MODEL  GEOM  -  1.6"  OIA  x  2.16"  LQ  RT  CYL  (STEEL) 
MODE  LWT- 0.904  LBS 
1  CALIBER  -  1  NOSE  PLAT  DIA 


FIGURE  8  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 
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FIGURE  9  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 
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NSWC  TR  81-50 


SHOT  #1711 

MODEL  QEQM  -  RIGHT  CIRC  CVL 
1  CALIBER  -  1  NOSE  FLT  DIA 
Po  -  1  ATM 

-  '  --  -  EXPERIMENTAL 
- -  -  THEORETICAL 


FIGURE  ie  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 
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SHOT  #1438 

MODEL  QEOM  -  RIGHT  CIRC  CYL 
1  CALIBER  >  1  NOSE  FLAT  OIA 
P,  -  1  ATM 

- -  EXPERIMENTAL 


FIGURE  17  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 
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SHOT  #1434 

MODEL  QEOM  -  RIGHT  CIRC  CYL 
1  CALIBER  >  1  NOSE  FLAT  DIA 
Po  -  1  ATM 

- -EXPERIMENTAL 

- -  THEORETICAL 


0.010  0.020  0.030  0.040 


TIME  (SEC) 


FIGURE  18  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 


SHOT  #1474 

MODEL  GEOM  =  STEPPED  CYL 
1  CALIBER  =  1  NOSE  FLAT  DIA 


FIGURE  19  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 
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(suaanvo)  U3.i3wvia  axiavd 


FIGURE  20  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 


#1419 


FIGURE  22  CAVITY  DIAMETER  AT  VARIOUS  DEPTHS  VS  TIME 


TIME  (SEC) 

diameter  at  various  depths  vs  time 


=  0  Ibi/aq  ft 


FIGURE  24  CAVITY  DIAMETER  AT  20  CAL  DEPTH  VS  TIME.  FOR  VARIOUS  SURFACE  PRESSURES 
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TERMS 


Missile  cA-oss-aectional  drag  area 
Lamina  rim  surface  area 

Cavity  drag  coefficient  of  projectile  based  on  Aj. 

Diameter  of  lamina  opening;  diameter  of  cavity  at  a  given  depth; 
missile  nose-flat  diameter 

Gravitational  acceleration 

Projectile  mass 

Atmospheric  surface  pressure 

Radius  (and  "depth")  of  spherical  lamina 

First  critical  minimum  depth 

Second  critical  minimum  depth 

Projectile  depth 

Diameter  of  axis-symmetric  projectile  body  at  aone  of  separation 

Elapsed  time  for  missile  to  reach  depth  r,  measured  from  entry 

Kinetic  energy  of  fluid  in  shell 

Kinetic  energy  of  projectile 

Projectile  velocity 

Projectile  velocity  at  water  impact 

Total  potential  energy  of  lamina/surface-pressure  system 
Potential  energy  of  lamina  due  to  gravity 


NSWC  TR  81'-59 

TERMS  (Cont'd) 


V£n  Initial  potential  energy  of  lamina/surface-preaaure  system  in  the 

at-rest  (i.e.,  6^  ■  0)  position 

Vp^aR  Potential  energy  of  lamina/surface-pressure  system  when  lamina  is  at 

peak  of  its  awing  (i.e.,  0^  ■  9„) 

Vpg  Potential  resulting  from  atmospheric  surface  pressure 

Wp  Work  done  by  pressure  Pq  acting  on  lamina 

0 

X,  Z  Space  saving  substitution  parameters 

a  Retardation  coefficient,  a  C[)A^P^/2mp 

^a*  Fluid  shell  displacement  angles 

9.  Initial  shell  angle 

0 

Maximum  possible  value  of  the  shell  angle  6^ 

Py  Density  of  water 
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uPPENDIX  A 

PROGRAM  LISTING  AND  SAMPLE  OUTPUT 


A-1 
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73/^4  OPT-1 


FTN  4. 0*452 


1 


5 


10 


15 


20 


25 


30 


35 


40 


45 


so 


55 


C 

C 

c 

c 

c 

c; 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

i; 

c 

t 

c 

t 

c 

c 

c 

c 

c 

L 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


-P»06kAm  OESCRIPTTOM 


THE  HHORRAM  LISTED  BELOW  IS  DESIGNED  TO  CASHY  OUT  A 
NUMERICAL  SOLUTION  OF  THE  THEOHETICAL  EQUATIONS  OF  MOTION 
WHICH  approximately  PBEDICT  BEHAVIOR  OF  WATER  ENTRY  CAV¬ 
ITIES  fOR  VERTICAL  WATER  ENTRY. 

THE  t'lUMlONb  ARE  TAKEN  FHOM  A  HYDRAULIC  CAVITY 
Flow  MOOEL*PROPOStU  by  G.  BIHKHOFF  and  R.  ISAACS  ("TRAN¬ 
SIENT  CAVITIES  IN  AIr  water  ENTRY,"  NAVORD  14V0  ,  1951), 
which  is  based  on  the  APPROXIMATE  HYPOTHESIS  THAT  STREAM¬ 
LINES  LIE  ON  CONCENTRIC  SPHERICAL  SURFACES  (OR  "LAMINAS") 

centered  at  the  point  of  impact. any  nose  configuration  for 

WHICH  THE  DRAG  COtFF.  IS  KNOWN  MAY  Rt  MODELLED. THF  EUUA- 
TIONS  AS  USED  IN  THE  PROGRAM  HAVE  BEEN  MODIFIED  TO  INCLUDE 
AN  ATMOSPHERIC  pressure  TERM. 

specifically  THE  PROGRAM  GENERATES  "(AMINA  PERFOR¬ 
MANCE  DATA  TABLES"  FOR  THE  LAMINAS  AT  USER  SPECIFIED  DE»THS 
(MEASURED  ALONG  ThE  LINE-OF-UHE)  ,BY  MEANS  OF  A  STEP-BY- 
STEP  lNTE(iRATION  OF  THE  LAMINA  MOTION  EONS . COMfl INING  THE 
UATA  in  THE  tables  FOR  LAMINAS  AT  SEVERAL  DEPTHS  MAKES  IT 
POSSIBLE  TO  OBTAIN  THE  PERFORMANCE  HISTORY  OF  THE  OVERALL 
CAVITY, THE  OUTPUT  UATA  ARE  DISPLAYED  IN  TFRMS  OF  POLAR 
COORoINATE  HOWEVER  FDR  CONVENIENCE  AN  APPROXIMATE  CAVITY 
diameter  is  COMPUItU  AND  DISPLAYED  HUT  IS  ACCURATE  ONLY 
FOR  sufficiently  small  values  of  THE  VARIAHLt  "ANGLE". 

DIRECTIONS  TO  SET  UP  THE  INPUT  DECK  APPEAR  BELOW 
ALONG  WITH  A  DESCRIPTION  OF  IMPORTANT  VARIABLES. THE  PRO¬ 
GRAM  FEATURES  A  PRINT  FREQUENCY  PARAMETER . "NSTORE" , WH I CH 

AIDS  IN  Reducing  (he  quantity  of  output  produced. also 

THE  USER  NEED  NOT  bUTHER  WITH  "MAX  ITERATION"  OR  "MAX  TIME" 
parameters  SINCE, when  THE  MOTION  OF  A  PARTICULAR  LAMINA 
IS  COMPLertO.COMPOIATIONS  FOR  HAT  LAMINA  ARE  TERMINATED 
AUTOMATICALLY. 

PROGRAM  AUTHUk-  M.  A.  METZGER  »  SEPT.  IVHfl 


4  input  OATA  * 


FIRST  CARD  - 

variables  ■  HEFDIA.RLNGTH.wEIGHT ,MMnRHO,DGAREA,UGCOEF 
DESCRIPTION  '  REFERENCE  D I  AM.  ( I N , )  ,  MOOFL  LENQTHdN.l, 

model  weight (lbs.) , MOD.  PENS.  (LB/CD. IN, )  , 
DRAG  AR£A(SU.IN.) fDRAO  COEFF..RESP. 

FORMAT  »  6F10.5 


SECOND  CARD 


VARIABLES 

DESCRIPTION 

FORlxt 


■  ENTVELtNOEPTH 

■  ENTRY  VELOCITY(FPS) ,  NUMHER  OF  SPECI* 
FIED  DEPTHS  PARAMETER,  PFSPECT 1 VELY . 

■  HO. 5, IS 


A-2 
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Ti/lh  0PT“1 


F7N  *.6*4ii2 


60 


6S 


7() 


7b 


eu 


ab 


90 


c 

t 

L 

c 

L 

c 

c 

L 

C 

L 

L. 

L 

C 

1. 

L 

c; 

«- 

c 

c 

c 

L 

C 

L 

C 

C 

c 

c 

(- 

L 

C 

I, 

c 

i; 

(- 

t. 


THIHl)  St  I  - 

VARlABLtS  ■  UtPTHtDELTAfNSTORE 
OtSCHlPTION  ■  tACH  CARD  THIS  SET  CONTAINS  A  USER 
iPtCIEIEO  DEPTH  A1.0N(3  WITH  VALUES  FOR 
THE  DESIRED  time  INCREMENT  AND  PRINT 
t'HEUUENCY  FOR  THAT  DEPTH. DEPTHS  SHOULD 
bE  IN  UNITS  OF  CALIdERS  RASED  ON  REt'OlA. 
SUBMIT  DNt  UF  THESE  CARDS  PER  SPECIFIED 
DEPTH. 

HiRmat  “  riFlO. 5*1*1 


«  VARIABLES  • 


VARIABLES  ■  ANtiLE  tANGVEL.  ANOACC 
DESCRIPTION  »  ANGLE  REPRESENTS  THE  ANGLE  SUHTENOtO  AT 
IHE  POINT  OF  entry  BY  THF  LlNE-OF-FIHE 
AND  THE  CAVITY  WALL  AT  DEPTH,  ANGVEL 
AND  ANOACC  ARE  ANGULAR  VELOCITY  ANU 
ACCELERATION  ASSOCIATTED  WITH  ANGLE. 
TYPES  ■  REAL 


VARIABLES  ■  WATRHO.pRESSH 

DESCRIPTION  "  WATER  DENSITY  AND  ATMOSPHERIC  PRESSURE 
TLRM  RESP. 
lYPES  =  REAL 


9b  L 

C 
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PROORAM  CAVITY  T.3//4  OPT«l 


FTN  A.6^452 


PROQRAM  CAVITY  1 1 NHUT » OUTPUT i TAPES- INPUT  * TAPfcb-OUTHUT ) 

C 

C 

t  ###»*#****#«##******##********#»*(i**#*#**#*****»**» 

c  preliminary  operations 

c  *##******#***«•***#**#****♦♦****♦*#**»•»♦*#•*»#*##* 

c 

REAL  MOMASS  .  MDLUEP 

dimension  depth  (lU)  «OELTA(10>  fNSTOREdO) 

QRAVTY»32.17ft 

ttATRHU-1.939A5 

PRESSR-211'S.H 

READ(S«100U)  REPDiAiRLNOTHtMElUHTfRMORHOiOOARcAiDGCOLI-' 
R£AD(5.101U)  ENTVtLiNDEPTH 
00  1  I-ltNDEPTH 

REAO<5»IOKO)  DEPTH(I) ,DELTA(I)*NST0RE(1) 

1  CONTINUE 
C 
C 

C  *  print  REFEhtNCE  INFO  TAdLE  *» 

C 

write (AtZOUO)  REFUIAtRLNOTHiOOCOEFtOOAREAtWEIGHTt 
.  RMOHHOtWAT«HO»PRESS8»8HAVTY«t'NTVfLiNDEPTH 

REFOlA-REFDIA/lS.O 
OOAHEA-DOAREA/HA.U 
M0MA9S-WEI0HT/32.I74 

RETARD- <DGLOeF*DGAHtA*WATRHO)/(2.0*MOv(ASS) 

WMITtf  (6,200&)  RETARD 

white (6*2010)  (OEPTH(I) ,I«1 *NDEPTH) 

DCONV-leO.U/3.1AIb92 

C 

c 

C  *  PRINT  TABLt  HEADIN5  * 

C 

WRITP  (6*20<e01 
L 
C 

c  ***##»*#--»****♦«*■**#*»***#*#*»»***«•»*##»****#**•* 

C  (iENL»^ATION  UF  LAMINA 

C  PERFORMANCE  DATA  TA9LES 

C  •**««•»•« *«««««*««*»**»*4>  »»««««»•*««*««••«««» »*«*»« 

c 

DO  So  I«1»NDFPTH 

c 

c 

C  •  INIIIALIZE  PARAMETERS  * 

C 

ATEMP-RETARO*MDMASS*ENTVEL*ENTVEL*EXP(-P,0*Rf:TAHD*DEPTH(  D* 
.  REFDIA) 

BTEMP-2, 0*3. 141692*  (DEPTH  <n**2)»(REFDlA**2) 

CTEMP-OEPTH  < I ) pREFOI A*WATRH0*GHAVT Y  +  PHFSSR 
COSMAX-1 . 0-ATEmp/ <BTEMP*CTEM») 

DEOANO>ACOS(COSMAX) *180.0/3.141692 

AHTlME“(tXP(RErAHO*l)EPTH<I )*REFDIA)-1.0)/(RFTARD*tNTVEL> 
WRITE (6* 20  30  I  UtPTHI I ) lOEL TA < I ) .NSTORE ( I ) tOERANG 
L 
L 


B0/1> 


A4 


C'  n  c:  Cl  cz  c:  cz  c:  a  cz  a  r.  r>  r:  n  r:  n  c'  c:  o  C"  c>  n-  ez  r  d  n  n  d  cz 
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PHO'.jHAM  CAVITY  73/74  0PT«l 


KTN  4.6*452 


*  commute  initial  AN3LE  * 

ANULE»ATAN(1.0/<2.0*0EOTHt  n ) ) 

CnSANB*COS(ANGLL) 

5lNANQ>5lN(AN0Lk.) 


“  CQNfUTE  INI  UAL  AN3.  VEL.  * 

Call  VhLOr  (PHEi>SH*WATRMO»ertAVTY,COSMAX.  ANGLF  . 
OEPlri  (  I  )  .WEFOl  A*AViGVF.L» 


*  COMPUTE  INI  UAL  AN3.  ACC.  • 

Call  ACCEL  (PWE!5bHiWATH-TO*URAWTY»COSMAX.AN3LF  * 

,  OEPTri(I) .REFDIA.ANGACC) 

MnLUEP  =  AL()8(HEI  ARD*ENTVEL*ABTIME  +  1  .01  /  (PETARn*REFDI  A) 
C  1  VUI  A*t>,  n*OEPTri  ( 1 1  *SI  NTANULE) 


4  HAOIAM  TO  Ufc.GREF  CONI/EHSION  * 

DFbANQ  =  ANr,LE*DCUNV 
UFOVEL ■ ANGVELAJCUNV 
0K(3ACC'AN(iACC«UCUNV 


4  PWliO  NEXT  LINE  OK  TAbUt  * 

W>^ITF.  (h. 20401  AbTIMF.iCAVDIA*OEGANU. 

.  UtOVEL*OEGACC,MULUEP 

STEP  BY  ultP  tNTEQHATION  OF 
MOTION  EOiNS  for  lamina  AT 
(.URRENT  OtPTHd) 

«<»•*•#«*** »*»*•«««» *«•»*«««»»«««»«»»«*•«««•»«•»«*« 

OLMANG=ANnLE 
I ST0«E=0 

ii;  lSTORE«lSTORE*i 


*  COMPUTE  NEW  ANGLE  ANO  ANGVEu  » 
AnGLE»OLOA.NG+ANGVEL«OEl.TA  (1 )  ♦0,5*ANOACC:*  (DELTA!  1 1**2) 

*  CHECK  for  NtO.  ANGuEflF  SO  « 

PROCEED  TO  NEXT  DEPTH 

lE-(ANQLt)  50*40.40 
•to  CONTINUt 

Call  VELOC (PHESSR,bATRRO*GRAVTY,COSMAX. angle. 

.  0EPTm(I),REF0IA,AN8VEL) 

IF  (ANGLE. LE.OLi/ANG)  ANOVEL—ANGVEL 


bO/1 


A-B 


n  r:  n  a  Ci  r  rirsr.  n 
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PROCHaM  C/ks/ITY  T3/74  OPT«l  FTN  4.6*452 


OLUANO-aNOLE 


*  Compute  new  ang.  acc.  * 

CALL  ACCEL (PREShH.WATRHOiORAVTYtCOS  4AX«AMOLF  « 
.  DLPTH(  I)  «REF0IA«A'46ACC) 

ABTIME-AHTIME*0ELTA( I ) 

IH  I  ISTOHE-t^STOREU  )  )  30.20.20 
■)il  CONTINUE 
STOP 


format  statements 


innu  format (6F10.S) 
lUlU  format (Flo. 5, 1B» 
lOi-U  format  (2Fld.S.  IS) 

2000  FORMAT ( IHl .///////.4bX. 

.  44M*  reference  information*. 

.  //////. ATX. 33HMO0EL  CONFIGURATION  «  SLIHROC  .//. 

.  44X.2AHM00EL  REF.  OIAM.  (IN.)  •  .Fa . 2 . / / . 4HX i 

.  22hmqi)el  Length  (in.)  «  .Fs.a./'/.sirx. 

.  20HCAV.  drag  COEFF.  «  . F5. 3. // ♦ A 3X . 

.  27MEFF.  OHaG  AREA  (SO. IN.)  ■  . K R .5 . // . A PX * 

.  22HMOnEL  WEIGHT  (LBS)  ■  .F8.5.//»AnX. 

.  aoHMOnEi.  density  (lhs/cu.in.)  ■  »ft.a.//.ai x. 

.  29H4ATEH  OtNS,  ( SL JG/CU. FT . )  ■  . FT , 5 . // , A2X  . 

.  26MATM.  PRESS.  (LOS/SQ.FT.)  ■  . F 7 . ? ♦ // . A I  X  . 

.  29HftHAV.  CONST,  (FT/3EC/SEC)  ■  ,  FF .  3 » /’/' ,  A  7  X  . 

.  23HENTHY  VtL.  (FT/SEC)  »  »F5 , 1 . 2/ .P6X . 

.  lAHNO,  depths  *  .12./) 

200S  F0RMaT(1H  .si  X.  IH-iHETAHD  (1/FT.)  ■  .Fb.6./I 

2010  format (IH  .AMX,21HDEPTHS  (CALIWERS)  «  . F A , ) . AFF . 1 . / ) 

2l)t  0  FORMAT  (IHl  .///////. 3bX, 

.bFH-  SPHERICAL  LAMINA  PERFORMANCE  DA 
.A  -»///', 53A.P9HZEHU  riME  REF.  »  WATER  IMPACT) 

2010  format (IH  .////////. 3AX.27H*  ***•*•**»#**#  , 

.  FA.1.38H  CAL  LAMINA  *****••♦#•**#  *,///. 

.  51«,19HTIME  STEP  (SEC)  *  .F6. 5 » / . FO X . 1 OHNSTOHE  »  , 

.  I3./.SIX.19HMAX  ANGLE  (DEG)  »  ,FS ,?.///-,  3TX  .AHAPS .  . 

.  8X.7HAPPROX,i iSX.AHANG. .RX.AHANG,  .PX.SHM0DEL./.37X, 

.  AHTJME.8X,bMCAy.0IA..2X,5HANGLE.7X,AHVEL..BX.AHACC.i 

.  8X.bHDEPTH./.36X»6H(SEC) »6X,7H (CAL"S ) . 3X ,5H ( DE(i ) . 

.  AX.9H (OEG/SLC) ,3X. 10H(DE3/SEC2) . AX,7H (CAU"S )  i/) 

20AO  FORMATdH  .3FX,F5.A.  8X  «F5 ,2  .  AX  .F5 ,2.  AX  ?  Fp  ,  2 , 1  X  ,  F  I A  .  1 . 3X  .  F7 . 2 ) 
END 


A-6 


BO/ll 
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SuBNOUri'iE  ACO.L 


ra/T't  oPTai 
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SUBHnUTINE  ACCEL  (HWES5H.«ATRH0..GHAVTY»C0SM AX  fANfiLE* 

.  UtPTHiHEFOIAtANGACC) 

COSANO'COSIANGLE) 

S1NAmS»SIni ANGLE) 

A*  l2.U-COSANO)*(i.O*COSAN3)/(COSAN()*<1.0-COSANfi)  ) 
HLOOX«AlOG(X) 

A»2.i)-4.0»COSANa 

B»RLf)liX*C0'3AMQ*  (2.U»COSAN3**3-2.0*COSANQ*r,0SANG-CO5ANU) 
Ca2.n«C0SA  ja/(SlNANQ*SINANG) 

Z= (A/ri) *0 

dummy 3  (  1  ,0»  (COSANo-COSMAX)  *Z)  /’(SINANe“WLOGX) 
ANGArC«-Z.U* (OERTrt*RtFOtA*GRAYTy*PRESSH/»<ATRHO) »DUMMYX 
.  ( (n£PTH*«EFUIA)**2) 

return 
END 


SURROUTIrE  Vkuuc 


73/74  OPT=l 


FTN  4.6*452 


5U9RnUTINE  VFLOC ( RhESSH . WATRMO . GRAY TY » CDSMAX • ANGLE ♦ 

.  UtPTH.REFDlA.ANGVEL) 

C0SA*i(5"C05  <  ANGLE) 

SINAmGbSIN (ANGLE) 

X» (2, U-COSANG) » ( 1 .U+COSANG) / (COSANG* ( I .O-COSANG) ) 
RLUGrALOG( A) 

ANGVFL  =  S0  >  I  (4.0*  (l'tPTH*HEFOl  A*Ct-)AVTY  +  OHE5SR/WATf(HO) 

.  It  (COSANG-COSMAX)  /  (  (  OEPTri*REFD  I  A  )  **?*S  I N  ANG 

.  i‘SINAI'tG*HLQG)  ) 

WLTUftN 

END 


B0/I( 


«0/l( 
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•REFERENCE  INFORMATION 

model  conk leUHATION  a  HIOHT  CYL. 

MODEL  ref.  01  am,  (IN->  •  i.sn 

MODEL  LENGTH  (IN.»  »  ?.l!i 

CAV.  DRAG  COEFF,  <*  .AOT 
tFF,  drag  area  (SO. in.)  »  1.76715 

MODEL  WEIGHT  (LOS)  »  1.07522 

MODEL  density  <LBS/CU,IN,)  »  ,2830 

WATER  OFNS.  ( SLUG/CU .FT  . )  >»  1,93945 

AIM,  PRESS,  (LBS/SU.FT.)  *  2116. HO 

GRAY.  CONST.  (FT/SEC/SEC)  *  32.174 

ENTRY  VEL.  (FT/SEC)  »  250.0 

NO.  DEPTHS  s  3 
retard  (1/FT.)  s  .28737 
DEPTHS  (CALIBERS)  -  10,0  20.0  30,0 


A-B 


NSA'C  TR  ft1~69 

-SPHEHICAL  LAMINA  PEPFOMmANCE  DATA 
ZEPO  TIME  REF.  «  WATER  IMPACT 


*»««••»•««*  «•»  10.0  CAL  lamina 


TIME  STEP  (sen  =  .00010 

NSTORF  a  S 
MAA  ANQLE  (DEG)  *  R.aJ? 


AHS. 

APPROX. 

ANG. 

AN6  . 

MODEL 

time 

CAV.OIA. 

ANGLE 

VFL. 

ACC. 

DEPTH 

(SEC) 

(CAL"S) 

(DEG) 

(DEG/SEC) 

(DE6/SEC2) 

<C6L''S) 

.nONO 

1.00 

2.86 

2524.23 

-2148373.6 

10.00 

.  0  0  g5 

37 

3.93 

1837.22 

-911064,9 

10.69 

.0070 

1.66 

4.75 

1486.39 

-548861,7 

11.36 

.0  0-75 

1.89 

5.43 

1257.18 

-386432.5 

12.02 

.0080 

2.10 

6.02 

1088.26 

-297266.2 

12.66 

.0065 

2.27 

6.53 

954.43 

-242139.3 

13.29 

.0090 

2.93 

6,97 

843.15 

-205324.0 

13.90 

.0095 

2.57 

7.37 

747,34 

-179366.2 

14.50 

.0100 

2.69 

7,72 

662.66 

-160333.9 

15,09 

.0105 

2.80 

8.03 

586.25 

-145974.0 

15.67 

.0  110 

2.89 

8.31 

516.15 

-134912.2 

16.23 

.0115 

2.97 

8.55 

450.94 

-126269.1 

16.78 

.0120 

3.05 

a.  76 

389.58 

-119469.8 

17,32 

.0  125 

3.11 

8.94 

331.25 

-114083.9 

17,86 

.0  130 

3.16 

9,09 

275.31 

-109861,4 

18.36 

.0135  ' 

3.20 

9.22 

221 .24 

-106594,0 

IB. 89 

.01  AO 

3.24 

9,31 

168.59 

-104140.8 

19.39 

.01A5 

3.26 

9,39 

116.96 

-102403.2 

19.89 

.illSO 

3.28 

9,43 

66.08 

-101314.7 

20.37 

.0155 

3.28 

9.45 

15.57 

-100835.6 

20,85 

.0160 

3.28 

9.45 

-28.66 

-100803.8 

21.32 

.0  165 

3.28 

9,42 

-76.22 

-101482.6 

21.78 

.0170 

3.26 

9,37 

-127.24 

-102697.1 

22.23 

.0175 

3.23 

9.30 

-179,02 

-10*571.7 

22.68 

.0180 

3.20 

9.19 

-231,93 

-107178.7 

23.12 

.0185 

3,15 

9.07 

-286.34 

-110624.3 

23.55 

.0190 

3.10 

8.91 

-342.71 

-115059.4 

23.98 

.0195 

3.03 

8.72 

-401.59 

-120696.5 

24.40 

.0200 

2.9b 

8.51 

-463.65 

-127035.9 

24.81 

.0205 

2.87 

8.26 

-529.74 

-136908.7 

25.22 

.0210 

2.77 

7.96 

-600.98 

-148547. B 

25.62 

.0215 

2.66 

7.66 

-678.86 

-163711.7 

26.01 

.0220 

2.54 

7.30 

-769.51 

-183911.0 

26.40 

.0  225 

2.46 

6.69 

-864.00 

-211650.6 

26.  79 

.0230, 

2.24 

6.43 

-979.11 

-251366.2 

27.17 

.0  235 

2.06 

5.91 

-1 118,70 

-311619.8 

27.54 

.0240 

1.85 

5.30 

-1297.04 

-411041.0 

27.91 

.0245 

1  .60 

4.60 

-1543,69 

-597968.7 

28.28 

.0250 

1.30 

3,74 

-1934.54 

-1040086.3 

28.64 

.0255 

.91 

2,60 

-2769.24 

-2802135.3 

28.99 

.0260 

.12 

.34 

-862860962.9 

29.34 
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»  •  #  #  ♦ 

«««»»• 

*»*  80.0  cal  lamina 

•o*  «****• 

*  «  •  «  ( 

TIME 

STEP  (SEC)  « 

.00010 

NSTORE  ■ 

5 

MAX  ANGLE  (t)EG1  " 

3.24 

AR5. 

APPROX. 

AN6. 

ANG. 

MODEL 

TIME 

CAV.OIA 

.  ANGLE 

VFL. 

ACC. 

depth 

(SEC) 

(CAL"S) 

(DEG) 

(DEGySECl 

(DEGX5EC2) 

(CAL"S) 

•  0U6 

1.00 

1.43 

760,50 

-455800.8 

20.00 

.nisi 

1.83 

1  .T7 

590.90 

-2S7576.3 

20.48 

.nis6 

1.42 

2.03 

434.92 

-176413.0 

20.96 

.nisi 

I.S7 

2.26 

408.34 

-133996.1 

21.43 

.ni66 

1.71 

2.44 

348.19 

-108615.4 

21.83 

.nirl 

1,82 

2.60 

298,29 

-92034.1 

22.34 

.nlT6 

1.91 

2.74 

855.27 

-80692.6 

22.78 

.01(^1 

2.00 

2,66 

217.08 

-72541.1 

?3.22 

•  Ola* 

8.07 

2.96 

13?. 38 

^66566 • 6 

23.65 

.0191 

8.12 

3.04. 

150.26 

-62136.2 

24.07 

.ni'96 

8.17 

3.11 

120.06 

-58354.9 

24.49 

.  nzol 

^  1 

3.16 

91.26 

-56469.1 

24.90 

.0806 

2.83 

3.20 

63,47 

-54B14.1 

25.31 

•  n8l  1 

8.8^) 

3.23 

36.35 

-53735.4 

25, T1 

.0816 

2.26 

3.24 

9.60 

-53321.9 

26.10 

.0881 

8.86 

3.24 

-14.93 

-53371.1 

26.49 

.0826 

8.85 

3.88 

-41,73 

-53944.5 

86.88 

•  0231 

8.83 

3.20 

-68.97 

-55092.4 

87.25 

.0836 

8.20 

3. IS 

-96.93 

-56833.6 

87.63 

.0841 

8.16 

3.10 

-125.97 

-59433.1 

88.00 

.0846 

8.11 

3.03 

-156.51 

-62920.4 

88.36 

.0861 

2. OS 

2.94 

-189,09 

-67622.3 

80.72 

.0  256 

1.96 

2.84 

-224.40 

-73971.5 

29.07 

.0861 

1.90 

2.78 

-263.44 

-38667.5 

89,4? 

.0266 

1.80 

2.57 

-307.64 

-94896.9 

29.77 

.oari 

1.66 

2.41 

-359.26 

-112815.1 

30.11 

.n8?6 

1,54 

2.21 

-42?. 07 

-140719.8 

30.45 

.0881 

1.36 

1.9B 

-503.15 

-188392.3 

30.70 

.0886 

1.19 

1.71 

-617,90 

-203058.5 

31.11 

.0891 

.94 

1  .35 

-809.73 

-532201.3 

31.43 

.0846 

.60 

.85 

-1298,35 

-1912956,5 

31.75 
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30.0  CAL  lamina 


time  step  (SEC)  ■ 
NSTORE  » 
MAX  angle  (DEG)  » 

.00010 

5 

1.48 

ARS. 

TIMF 

(SEC) 

APRROX. 

CAW.OIA.  ANGLE 
(CAL"S)  (DEG) 

ANG. 

Vf  L. 

(DEG/SEC) 

ANG. 

ACC. 

(DE6/SEC2) 

model 

DEPTH 

(CAL''S) 

0270 

1.00 

0275 

1.13 

0280 

1.24 

0286 

1.33 

0290 

1.40 

0295 

1.46 

0300 

1.50 

,0305 

1.53 

.0310 

1.55 

.0315 

1.55 

.0320 

1.55 

.0325 

1,53 

.0330 

1  .50 

,0335 

1.46 

.0340 

1.40 

.0345 

1.33 

.0350 

1.25 

.0365 

1.14 

.0360 

1.01 

.0  365 

.83 

.0370 

.60 

.0375 

,15 

.95 

28».e2 

1.08 

230.10 

1.19 

166.17 

1.27 

150.56 

1.34 

120.12 

1.39 

93.05 

1.43 

68.22 

1.46 

44.87 

1.48 

22.42 

1.48 

.41 

1.48 

-21.61 

1.46 

-44.03 

1.43 

-67,34 

1.39 

-92.10 

1.34 

-119.08 

1.27 

-149.37 

1.19 

-184.76 

1,09 

-228.29 

.96 

-286.29 

.80 

-374.45 

.57 

-562.06 

.14 

-2170.53 

-140130. A  30,00 

-•JVfcSA.S  30. 3A 

-78160.5  30.67 

-65303.5  31.00 

-57071.8  31.33 

-61608.7  31.65 

-47969.1  31.97 

-45637.5  32.29 

-44328.6  32.60 

-43897,3  32.91 

-44297.9  33.21 

-45572,7  33.51 

-47863.0  33.81 

-51447.5  34.11 

-56830.8  34.40 

-64936.6  34.69 

-77564.2  34.97 

-98613.0  35.26 

-138067.7  35,54 

-229697,3  35.82 

-675360,9  36.09 

-31522971.7  36,36 
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appendix  » 


nPRTVATION  OF  LAMINA  MOTION  EQUATIONS  WHICH  INCLUDE 

derivation  ot  surface  pressure 


LINE  OF  FIRE 


NSWC  TR  81-39 


Thii  derivation  uaee  an  energy-approach  to  obtain  the  equations  governing 
the  motion  of  the  laminai  Refer  to  the  figure  above  throughout  the  derivation. 

First,  obtain  expression  for  the  potential  due  to  surface  pressure  Pq, 
which  is  equal  to  negative  of  the  work  done  by  Pq  moving  thru  the  displacement 
(9b  -  90°): 


.  dVp  -  -  dWp  -  -  Pq  A,  r  dOj,)  (1) 

•^o  *^0 

where  the  area  Ag ,  as  shown  in  the  figure,  is  given  by: 

Ag  ■  2Tfr  Ar  sin  9b 


Substitute  expression  for  Ag  into  Equation  (1)  and  integrate  to  obtain  an 
equation  for  Vp  in  terms  of  8b,  thus: 


Pq  211  r2  Ar  sin  9b  d0b 
90° 


2it  r2  Ar  p^  cos  0  b 


90° 


or: 


Vp  »  -  2tT  r2  Ar  Pg  cos  0b  (2) 

0 

It  can  be  shown  that; 

cos  0  b  “  9  a  "  ^ 


Substituting  this  into  Equation  (2)  yields  the  expression  for  Vp  in  terms  of 
angle  ° 


Vp  -  -  2  r'^  r  Pg  (cos  g  -  1)  (3) 

Next  consider  gravity  potential  Vg.  If  the  original  undisturbed  water  surface 
is  taken  as  the  reference  datum  for  zero  gravity  potential  then 
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Vg  ■  -  2ti  r2  Ar  p^,  g  h 
where,  as  shown  in  figure: 

ii  “  r(2  COB  00  -  l)/2 

Hence  gravity  potential  in  terms  of  shell  angle  is 
Vg  "  TTr^  Py  Ar  g  (1  -  2  cos  O^) 

Finally,  combine  potentials  due  to  gravity  and  surface  pressure  to  obtain  the 
total  potential  of  the  shell/pressure  system 


V  -  Vg  4 


3 

TTr  p  Ar  g  (1  -  2  cos  Q  )  -  2TT  Ar  P  (cos  0  -  1) 

V  A  O  A 


(A) 


Next,  consider  tlie  kinetic  energy  of  the  moving  shell, 
shell  angle  0^,  the  expression  for  the  kinetic  energy  can  be 


In  terms  of  the 
shown  to  beA“i 


Tf 


^  r^  Py  Ar  00  sin^  6a  In 


(2  -  cos 


(1  +  cos 


V 


cos  0.  (1  -  cos 


9.^ 


(5) 


For  the  total  energy  in  the  conservative  shell/surface-pressurG  system  the 
following  relation  holds  at  all  times 


(Tf  +  V)i  -  (Tf  +  V)2  -  Const. 


(6) 


Wtien  the  shell  is  at  the  peak  of  its  swing,  its  kinetic  energy  is  zero  and  the 
total  potential  of  the  system  is  at  a  maximum  which,  from  Equation  (6),  is  equal 
to  the  total  energy  of  the  shell  at  any  other  time. 

(Tf  +  V)  -  V(peak)  .  (7) 


At  the  peak  position  6q  is  also  at  a  maximum  and  is  designated  as 
the  maximum  value  that  9^  ever  reaches.  When  the  previously  derived 
expressions  for  Tf  and  V  are  substituted  into  Equation  (7)  and  the  shell  angle 
in  the  expression  for  V(peak)  is  replaced  by  6,^,  Equation  (7)  becomes 


A'iAbelson,  H. ,  "Cavity  Shapes  at  Vertical  Water  Entry  -  A  Comparison  of 
Calculated  and  Observed  Shapes,"  NOLTR  67-31  (1967). 
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Y  f  ^  Pw  [^]  ■•■  11  (Ar)  p4,  g  (1  -  2  coi  0^) 

~  211  r2  (A:)  P©  (cos  0^  “  1) 

■  H  r^  py  CAr)  g  (1  -  2  coi  0^)  -  2ti  (Ar)  P©  (coi  0„  -  1) 


where 


X  »  (2  -  COB  6j)  (1  +  COB  6||)/coa  0^  (1  -  coi  0^) 


Solving  this  lait  expression  for  the  angular  veloctiy  0©  yields  the  governing 
differential  equation  for  the  shell  angle 


5. 


n 

COB  6  -  COB  6 

a  m 

iin^  0  In 

s 

’(2  -  cos  Oa)  (1  +  COB  0a)n 

cos  0a  “  cos  0a)  J 

(8) 


The  expression  for  angular  acceleration  0^  can  be  found  by  differentiating 
Equation  (8)  with  respect  to  time 

-2  |rg  (Pg/Pyj]  f  ^  0,  ”  ‘^0®  tZ)"] 

8.  ■  ■  -J— ?! -  [  nreTulVl - 


(2  -  cos  6a)  (1  cos  Og) 
COB  9a  (1  “  COB  9a) 


(9) 


Z  ■ 


2  -  if  cos  9a 


2  cos  0a 


cos  9^  (cos 3  0^  -  2  cos^  8,  -  cos  0^  +  2)  In  [xj  gi©'/ 


The  expression  for  cos  0^  is  derived  next.  At  the  peak  of  its  outward  swing 
(corresponding  to  0,  “  9^)  all  the  kinetic  energy  imparted  initially  to  the 
lamina  by  the  passing  projectile  has  been  converted  entirely  into  potential 
energy  thus: 

'^peak  ”  '^in  “ 

where  Tp,  the  kinetic  energy  of  the  projectile »  is  given  by 


^  ^  2 
Tp  “  Y  "'p 
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If  Che  cei/iCy  drag  coefficient  ia  asaunmad  Co  be  conaCanC  then  Tp  can  be 
written  in  terms  of  lamina  depth  r  at 


r  nip  Vp2  e 
^  0 


-2ar 


The  expreationa  for  Tp,  Vp^^l^  and  are  substituted  into  Equation  (10) 
and  Che  differentiation  carried  out  to  obtain 


It  p  {li)  g  (1  -  2  cos  0n,)  -  (-it  r^  p  (Ar)  g) 


w 


“  211  (Ar)  Pq  (cos  9m  “  1)  ■  ot  mp  Vn  e  ^°*’^(Ar) 


P  "Pr 


Solving  for  cos  6„,  yields  the  final  desired  relation 


cos  9m  “  2 


am  Vn  e 
__J>L.o 


-2ar 


IT  r2  (r  g  +  P^)_| 
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